Unusual hyperfine interaction of Dirac electrons and NMR spectroscopy in graphene 
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Theory of nuclear magnetic resonance (NMR) in graphene is presented. The canonical form of 
the electron-nucleus hyperfine interaction is strongly modified by the linear electronic dispersion. 
The NMR shift and spin-lattice relaxation time are calculated as function of temperature, chemi- 
cal potential, and magnetic field and three distinct regimes are identified: Fermi-, Dirac-gas, and 
extreme quantum limit behaviors. A critical spectrometer assessment shows that NMR is within 
reach for fully ^^C enriched graphene of reasonable size. 
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Nuclear mametic resonance (NMR) is a powerful spec- 
troscopic tool [l| and an architecture for quantum infor- 
mation processing [3, Q at the same time. Both of these 
applications arc possible due to the relatively weak in- 
teraction of the nucleus with its environment. This weak 
interaction is sufficient to probe the electronic state of its 
vicinity, which yields information about the local elec- 
tron bonds or about the correlated behavior of electrons 
as e.g. in superconductors NMR quantum comput- 
ing exploits that the nuclei are well isolated from the 
environment thus there is a longer time window for the 
manipulation and detection of the nuclear quantum state. 

For both kinds of applications, the important NMR 
parameters are the shift of the NMR resonance with re- 
spect to a standard, and the decay of the longitudinal 
magnetization to its equilibrium value, the spin-lattice 
relaxation time, Ti. These were extensively studied in 
solid state systems both theoretically and experimentally 
[J, However, the body of NMR experiments were fo- 
cused on three-dimensional systems which stemmed from 
the unavailability of stable, inherently two-dimensional 
materials. The discovery of graphene, a single stable 
sheet of carbon atoms in a hexagonal lattice [6| , enables 
studies of an exactly two-dimensional system. Its quasi- 
particles follow a linear band dispersion, causing the elec- 
trons to behave as massless Dirac fermions, which gives 
rise to unique transport and magnetic properties Sim- 
ilarly, unusual electron-nuclear interaction is expected. 

In a metal, the NMR measurables are most affected by 
the surrounding electrons through the electron-nuclear 
hyperfine interaction (HFI). The standard, text-book 
form of the HFI of nuclei and conduction electrons leads 
to the Hamiltonian i/nFi = ^^orb + ^^spin 01 : 
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magnetism, the second (Hspin) contains the electron spin- 
dipole interaction and the so-called Fermi-contact inter- 
action. ^0 is the permeability of free space, 7„ is the 
nuclear gyromagnetic ratio, I is the nuclear spin, g ~ 2 
is the g- factor of the electrons. S, p, and r are the elec- 
tron spin, momentum, and vector operators, /ib is the 
Bohr magneton and /.j^ = m/m*^^ is the effective or- 
bital Bohr magneton y, where to* is the effective band 
mass and m is the mass of a free electron. 

At first, it is not obvious how to generalize the or- 
bital term of the HFI to massless Dirac fermions and its 
derivation is one of the primary goals of this work. Sec- 
ond, the unique properties of the conduction electrons of 
graphene are expected to give rise to unique relaxation 
and NMR shift behaviors. E.g. deviation from the Kor- 
ringa relation [l[ , that is an important benchmark of non 
Fermi-liquid behavior, is expected. 

Here, we show that the canonical description of the 
hyperfine interaction is modified for the massless Dirac 
fermions and we derive the hyperfine Hamiltonian paying 
special attention to obtain the appropriate orbital contri- 
bution. We identify different regimes based on the NMR 
measurables; Fermi-, Dirac-gas, and extreme quantum 
limit behaviors. We also discuss the feasibility of bulk 
NMR spectroscopy on graphene with a critical evalua- 
tion of NMR spectrometer performance. 

The low energy excitations in graphene are described 
by the two-dimensional Dirac equation 0] : 
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Here, the first term (i/orb) is due to the electron orbital 



where w 10*' m/s is the Fermi velocity of graphene, 
and the pseudospin variables (Pauli matrices, a) spring 
from the two-sublattice structure. The HFI in graphene 
is derived following Abragam Q by treating the nucleus 
as a magnetic dipole with m = fi-7„I. Its vector potential, 
A(r) = ^ "^j , is inserted into the kinetic momentum 
as p p + eA in addition to the electron and nuclear 
Zeeman terms. This calculation gives the effective HFI 
Hamiltonian in graphene as 
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where j = evpcr is the electric current operator in 
graphene and cr is a vector of the Pauh matrices. The 
first term describes the interaction of the nucleus with 
the orbital motion of the Dirac electrons 0, which con- 
tains Iz only as the electrons are confined in the plane. 
The spin-dipole and Fermi-contact terms are unchanged 
with respect to their usual forms. 

The orbital term in Eq. ([3]) differs significantly from 
the usual form in Eq. ([T]) as the orbital magnetic moment 
(r X j) replaces the usual term {^j^r x p). This is the 
result of the peculiar form of the current operator for 
Dirac electrons j ~ <t, which is also is responsible for the 
jittery motion of the center of mass coordinate known as 
Zitterbewegung 

Eq. (P) can be deduced formally 
from Eq. ^ by using j ~ ep/m* for a normal metal. 

A unique property of the orbital magnetic moment of 
graphene is that it remains invariant in an applied mag- 
netic or gauge field, since j is insensitive to the vector 
potential. We mention that the proper orbital angular 
momentum of Dirac particles is still r x p in the sense 
that it is responsible for rotations in the x ~ y plane, 
which differs from the orbital magnetization. We also 
note that there are no higher order terms in the vector 
potential in the graphene HFI Hamiltonian due to the 
linearity of the Dirac equation. 

The second quantized form of the orbital part of the 
interaction in graphene is obtained as 
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where Joib = /^o^7nefF/2Ac, and /(k, k', a, a') = 
{aa' - cxp[i{(pk - (pk')]){Q:k + a'k')/2\k-k'\, and c^^^ 
creates a quasi-particle with energy Ea(\i) and real spin 
s, ifk is the angle of k with the axis, is the unit 
cell area, and N is the number of unit cells. The interac- 
tion is bounded as |/(k, k', a, a')| < 1. The magnitude 
of the orbital term is estimated as Jorb ~ 21 MHz using 
7(i3C)/27r = 10.7 MHz/T. 

The effective interaction describing the hyperfine in- 
teraction in graphene is obtained from Eq. [31 as 
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where A is a 3x3 tensor with diagonal elements. Of these, 
the traceless ones are due to the spin-dipole interaction as 
Adip{x,y) : Adip{z) = —A^ip : 2A<iip and the scalar term, 
Aiso, is given by the isotropic Fermi-contact interaction. 



First principles calculations [llj gave A^ip = 73 MHz and 



Aiso = -44 MHz, which gives (-117,-117, 102) MHz for 
the diagonal elements of A. We note that the first prin- 
ciples value of Adip agrees well with the A^ip = 91 MHz 



obtained for the Pz orbital of a free carbon atom [l|, [l^] , 
which confirms that it is indeed the relevant orbital in 
graphene. 

Upon establishing the hyperfine interaction in 
graphene, we turn to the calculation of the NMR mea- 
surables. For a given magnetic field, terms of Eq. [S] 
perpendicular and parallel to the field contribute to re- 
laxation and to the Knight shift, respectively [l,[l3|. The 
spin-lattice relaxation rate, 1/Ti and the Knight shift, K, 
for a given magnetic field direction {i = x,y, z) are |14| 
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with Cf = + Si.,z'2Jlh)' 7e is the gyromag- 

netic ratio of electrons, p{E) is the quasi-particle density 
of states (DOS), and fi is the chemical potential. Ti and 
K for an arbitrary field direction is readily obtained by 
angular dependent combinations [isj . 

The orbital interaction involves only Iz , thus it affects 
Ti only when the field is in the graphene plane (z = 
X, y), which explains the 2J^^^ term (the factor 2 comes 
from the spin degeneracy). The orbital term does not 
contribute to the Knight shift even for a magnetic field 
along z in a manner analogous to demagnetization. The 
spin part of the HFl contributes to a '--^ 15 % anisotropy 
of Ti for in and out of plane magnetic fields but the 
orbital term makes it nearly isotropic. More accurate 
statements require the first principles calculation [llj of 
Jorb- The Knight shift changes sign and drops by 15 % 
from in plane to out of plane fields. We omit the i index 
from Ci in the following. 

We distinguish two scenarios for the DOS in the fol- 
lowing calculation: (i) absence of Landau levels and 
(ii) where the presence of Landau levels is important. 
Scenario (i) occurs for three cases: when magnetic 
field is in the plane, when magnetic field is arbitrary 
but level broadening due to F or T makes the Lan- 
dau levels undistinguishable around p (the criterion is 
VpeBjj, <max(r, fefiT)), or in the vicinity of the DP 
point (i.e. fx is small) when the lowest Landau level 
is significantly broadened due F or T (the criterion is 
vpVieBh <max(r, /cbT)). 

For scenario (i), the magnetic field-free DOS can be 
used in the calculation and it reads as: 
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FIG. 1: (Color online) The nuclear spin relaxation rate (main 
figure) and the Knight shift (inset) are shown with F = 0.1 -El 
and D = 1000 -El as a function of the chemical potential. The 
blue solid/red dashed line refers to the presence/absence of 
magnetic field at T = 0, the black dash-dotted line corre- 
sponds to ksT — -El in the presence of magnetic field. In- 
creasing fi OT T makes the Landau level structure disappear, 
and scenario (ii) is replaced by scenario (i). 



per spin and C atom. The resulting relaxation rate is 
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which increases as max(/-t^, ('Kk-QT)'^ / . Away from the 
Dirac point, dominates, and the temperature becomes 
important only near the DP. 

Exactly at the DP, Ti diverges as Ti ~ T^^, therefore 
the nuclear spins are not relaxed by conduction electrons 
at T = due to the absence of charge carriers at the 
charge neutrality point. In the presence of impurities, 
the DOS at the DP reads as [11 
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with F the scattering rate and D the cutoff in the con- 
tinuum theory. Therefore, the aforementioned divergence 
of the clean system weakens to Ti - (F^ ln^(£)/F)r)-\ 
reproducing the Fermi-gas behavior. Since the DOS is 
finite at the DP due to impurities, the Dirac nature of 
the quasi-particles is lost at this level. 
The Knight shift is evaluated as 
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It can be approximated by -fC ~ max(2A:BTln2, |/i|). Im- 
purities provide a finite DOS even at the DP, therefore 
the Knight shift stays finite there as if ~ Fln(-D/F). 

Concluding scenario (i), we give 1/TiT for the case 
of chemical doping of, or chemisorption on the graphene 
layer. E.g. for an AC^- composition, where A is an alkali 
atom with full charge transfer, there is an extra 2/xAc 
electron density to each lattice site. This translates to 



a chemical potential shift of /i 
leads to a relaxation rate as 



hvp^2'K / xAc, which 
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or Ti w 500 (sK) • x/T. This gives Ti « 10 s at 300 K 
for x=8, that is a usual doping level for graphite [lB|. It 
shows the sensitivity of the NMR properties for doping 
or chemisorption, which may lead to a sensor application 
of graphene. The chemical potential can be also tuned 
by gate voltage with a less dramatic effect on Ti. 

For scenario (ii). Landau level formation is important, 
and the continuous spectrum is replaced by discrete Lan- 
dau levels as Ena = ctEj^V^ where a = ±, n is non- 
negative integer, -El = vp\/2heB^ is the Landau scale, 
and Bz is the perpendicular component of the magnetic 
field. With this, the DOS reads asfH 
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where ^/{x) is Euler's digamma function, and reduces to 
Eq. ([7|) for clean systems with vanishing magnetic field. 
The l/TiT and the Knight shift are shown in Fig. [T] 
separately as a function of the chemical potential. These 
display the characteristic de Haas-van Alphen like oscil- 
latory behavior at high magnetic field and low p, and T, 
that we refer to as the extreme quantum limit (EQL). 

Calculation of the relaxation rate and Knight shift al- 
lows to test the validity of the Korringa relation, i.e. 
whether 1/TiTK^ = const, holds. In general, the Ko- 
rringa relation is valid for a Fermi-liquid. In particu- 
lar for a non-interacting Fermi-gas [22] {l/TiTK'^)-p — 
47rfcB(7n/7e)^/^- For graphene within scenario (i) and 
in the limit of {jjL^k-QT) ^ F, which is referred to as the 
scaling limit, it reads as 
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which depends only on the ratio of fi and T, and F{x) is 
a universal scaling function 
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which is even in x and satisfies F(0) = TT^/3ln{A) w 1.71 
and -F(oo) — 1, and is shown in Fig. [2l For in-plane field, 
C'^/A'^ < 1, while for perpendicular field C^/A"^ > 1. 
For /.J ^ k^T, the DOS is finite, and nothing distin- 
guishes graphene from a conventional metal because only 
one branch of the "V" -shaped dispersion is seen due to 
the smallness of T, therefore the usual Korringa relation 
is satisfied. In the opposite limit, (/i <C k^T), the Kor- 
ringa relation leads to a constant, -F'(O), times bigger than 
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FIG. 2: (Color online) The Korringa relation, normalized to 
its Fermi-gas value, as function of n/ksT (left panel) and 
fi/E-L (right panel). The usual Korringa relation is recovered 
for ksT but for increasing T the normalized Korringa re- 
lation increases and saturates to 7r^/31n^(4). The right panel 
shows the Korringa relation for F = 0.1 -El, which allows the 
visibility of the lowest Landau levels. The temperature is var- 
ied as ksT/Ei^—O.OS, 0.1, 1 and 20 from bottom to top. For 
ksT ~ F, separate peaks indicate the Landau level structure, 
and the curves cross over to the field-free scaling limit with 
increasing T. 



its conventional value, which signals the nature of Dirac 
fermions. The crossover can be explored even away from 
the DP by fixing the chemical potential to a finite value, 
and sweeping the temperature. Right at the DP, impu- 
rities spoil the crossover and re-establish the Fermi-gas 
relation for (/cbT, ^) <C P. 

The Korringa relation can be numerically evaluated 
in the presence of Landau levels (i.e. for scenario (ii)) 
using Eq. [T^] and Fig. O shows the result. For small 
T and F, the oscillatory behavior due to Landau levels 
in the DOS characterizes the Korringa relation. When 
(^bT, F) > El,, the Landau levels are smeared and the 
magnetic field does not play an important role thus the 
scaling limit is restored. 

In Fig. [3l we summarize our findings on the NMR 
properties in the form of a "phase diagram" . The extreme 
quantum limit shows up only at low temperatures and 
small chemical potential, when the Landau level struc- 
ture is visible. Larger F, i.e. presence of defects, favors 
the Fermi-gas region. 

We finally comment on the feasibility of NMR experi- 
ments in graphene. NMR is known to have a low signal 
sensitivity albeit its tremendous utility. In graphene, the 
NMR active ^'^C nuclei has a low abundance (c = 1.1 %) 
and a low gyromagnetic ratio, 7(^'^C) ~ 7(^H)/4. NMR 
spectrometers are characterized by the limit of detection 
(LOD) parameter, i.e. the number of nuclei required 
for a signal-to-noise ratio of three i n a single acquisi- 
State-of-the-art spectrometers 



tion. State-of-the-art spectrometers [17| have LODo = 
10^^/a/Hz for spins with sample and detector at 300 
K in a 14 T magnetic field (z^(^H) = 600 MHz). For a 
general case the LOD is [isf : 



FIG. 3: The schematic "phase diagram" of graphene accord- 
ing to NMR, the phases are separated by solid lines. The 
boundaries denote smooth crossovers. The solid lines move 
to the dashed ones for increasing disorder. 
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Here Tg* the apparent decay time of the NMR time- 
domain signal which contains the spin-spin relaxation 
time, T2, and the magnetic field inhomogeneity due to 
defects and the magnet. The Curie T dependence of the 
NMR signal is described by the sample temperature, T,. 
NFi-ci is the receiver noise factor relative to a receiver at 
300 K. 

Clearly, low sample temperature, low detector noise, 
and highly ^'^C enriched graphene are required for an 
NMR study. Sample temperature down to 1 K is cus- 
tomary in solid state NMR and NFj-ei = 1/8 was reported 
for cryo-probe NMR llQll . We estimate from NMR data 
on graphitic carbon [H [1^ that the FWHM = I/ttT^ 
is 50 ppm at 14 T of fully ^^C enriched graphene giv- 
ing T2 = 10 /iS. This estimate assumes either a sin- 
gle graphene sheet or a set of graphene layers oriented 
alike. These factors give an LOD for ^'^C graphene of 
8 • 10^^ which corresponds to a surface of 0.63 mm^. 
We think that synthesis of a fully ^^C isotope enriched 
graphene with such an area (not necessarily of a single 
piece) is within reach as fully ^^C enriched graphite was 
recently synthesized We expect that a dedicated 

NMR microcoil setup, prepared by lithographic methods 
[tH I would further decrease the LOD value and the re- 
quired graphene sheet area. 

In summary, we generalized the canonical theory of hy- 
perfine interaction between nucleus and conduction elec- 
trons for graphene. The orbital part of the HFI differs 
from its usual form as it does not involve the angular mo- 
mentum. We identified three distinct regimes in graphene 
based on the NMR measurables: Fermi- and Dirac-gas 
phases, and the extreme quantum limit. We argue that 
NMR on graphene is within realistic reach. 
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